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We consider eigenvalue problems 



Q . (1.1) — w" + P{z, ai)w = Xw, y{z) — )■ as 2; — )■ oo, z E LiVJ L2. 

^-^ ■ Here P is a polynomial in the independent variable z, which depends 

on a parameter a, and Li, L2 are two rays in the complex plane. The 
set of all pairs (a, A) such that A is an eigenvalue of (11. ip is called the 

r> I spectral locus. 

c^ • Such eigenvalue problems were considered for the first time in full 

generality by Sibuya [IQ] and Bakken [2]. Sibuya proved that under 
certain conditions on Li, L2 and the leading coefficient of P, there 
exists an infinite sequence of eigenvalues tending to infinity. If 

(1.2) P{z, a.) = z'^ + ad-iz'^~^ + . . . + aiz, 

where a = (ai, ..., a^.i), then the spectral locus, which is the set of 
all (a, A) G C^ such that A is an eigenvalue of ll.H is described by an 
equation P(a, A) = 0. Here P(a, A) is an entire function of d variables, 
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called the spectral determinant. So the spectral locus of f ll.ip . (11. 2p 
is an analytic hypersurface in C^. It is smooth [2] and connected for 
d>3 [HISQ]. 

In the first part of this paper we study what happens to the eigen- 
values and eigenfunctions when the leading coefficient of P tends to 
zero. 

Bender and Wu [8] studied the quartic oscillator as a perturbation 
of the harmonic oscillator: 

(1.3) -w" + {ez^ + z^)w = Xw, ti;(±oo) = 0. 

Here and in what follows w(±oo) = means that the boundary condi- 
tions are imposed on the positive and negative rays of the real line. It 
has been known for long time that the eigenvalues of (11.31) converge as 
e — ;■ 0-|- to the eigenvalues of the same problem with e = 0, but they 
are not analytic functions oi e a.t e = (perturbation series diverge). 
To investigate this phenomenon, Bender and Wu considered complex 
values of e and studied the analjd;ic continuation of the eigenvalues as 
functions of e in the complex plane. Their main findings can be stated 
as follows: the spectral locus of the problem (11.31) consists of exactly 
two connected components; for e ^ 0, the only singularities of eigen- 
values as functions of e are algebraic branch points. These statements 
were rigorously proved in [TH]. Discoveries of Bender and Wu gener- 
ated large literature in physics and mathematics. For a comprehensive 
exposition of the early rigorous results we refer to [H] . 

To perform analytic continuation of eigenvalues of (II. 3p and similar 
problems for complex parameters, one has to rotate the normalization 
rays where the boundary conditions are imposed. One of the early 
papers in the physics literature that emphasized this point was [7]. 
Thus physicists were led to problem (II. ip . previously studied only for 
its intrinsic mathematical interest. 

An interesting phenomenon was discovered by Bessis and Zinn- Justin. 
For the boundary value problem 

—w" + iz^w = Xw, w{±oo) = 0, 

they found by numerical computation that the spectrum is real. This 
is called the Bessis and Zinn- Justin conjecture (see, for example, his- 
torical remark in [5]). This conjecture was later proved by Dorey, Dun- 
ning and Tateo [151 ESj with a remarkable argument which they call 
the ODE-IM correspondence, see their survey [17]. Shin [37] extended 
this result to potentials 

(1.4) — w" + (iz^ + iaz)w = Xw, w{±oo) = 0, 
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with a > 0. 

These results and conjectures generated extensive research on the so- 
called PT-symmetric boundary value problems. PT-symmetry means 
a symmetry of the potential and of the boundary conditions with re- 
spect to the reflection in the imaginary line z i— >■ — ^. PT stands for 
"parity and time reversal" . 

It turns out that the spectral determinant of a PT-symmetric prob- 
lem is a real entire function of A, so the set of eigenvalues is invariant 
under complex conjugation. In contrast to Hermitian problems where 
the eigenvalues are always real, the eigenvalues of a PT-symmetric 
problem can be real for some values of parameters, but for other values 
of parameters some eigenvalues may be complex. So we can see the 
"level crossing" (collision of real eigenvalues) in real analytic families 
of PT-symmetric operators, the phenomenon which is impossible in 
the families of Hermitian differential operators with polynomial coeffi- 
cients. 

In this paper, we first consider the general problem fll.ip and the 
limit behavior of its eigenvalues and eigenfunctions when 

(1.5) P{z) = tz''\amZ^\v{z). 

with d > m > degp, as t — > 0, while the coefficients of p are restricted 
to a compact set and am does not approach zero. Then we apply 
our general results to certain families of PT-symmetric potentials of 
degrees 3 and 4, and prove some conjectures made by several authors 
on the basis of numerical evidence. 

In particular, our results for the PT-symmetric cubic 01.4p imply 
that no eigenvalue can be analytically continued along the negative 
a-axis, and the obstacle to this continuation is a branch point where 
eigenvalues collide. 

Another result is the correspondence between the natural ordering of 
real eigenvalues of (II. 4p for a > and the number of zeros of eigenfunc- 
tions that do not lie on the PT-symmetry axis, conjectured by Trinh 
in IH]. This correspondence is similar to that given by the Sturm- 
Liouville theory for Hermitian boundary value problems. 

A different approach to counting zeros of eigenfunctions is proposed 
in [27], where the authors prove that for a large enough, the n-th 
eigenfunction has n zeros in a certain explicitly described region in the 
complex plane. 

The plan of the paper is the following. In Section 2 we prove a 
general theorem on the continuity of discrete spectrum at t = for 
potentials of the form (II. 5p . with boundary conditions on two given 
rays. Previously such problems were studied using the perturbation 
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theory of linear operators in [261 El [II]- Our method is different, it is 
based on analytic theory of differential equations. 

Verification of conditions of our general result in Section 2 is non- 
trivial, and we dedicate the entire Section 3 to this. The question 
is reduced to the study of Stokes complexes of binomial potentials 
Q{z) = tz'^ + cz^^ d > m, which is a problem of independent interest, 
so we include more detail than it is necessary for our applications. The 
Stokes complex is the union of curves, starting at the zeros of Q, on 
which Q{z) dz^ < 0, so they are vertical trajectories of a quadratic dif- 
ferential. Stokes complexes occur in many questions about asymptotic 
behavior of solutions of equations (II. ip . Our study permits us to make 
conclusions on the behavior, as t — ?■ 0, of the Stokes complexes of po- 
tentials P{z) = tz''" + amit)z"^ + pt{z) where am(t) — )■ c 7^ and pt is a 
family of polynomials of degree m — 1 with bounded coefficients. We 
mention here [35] where a topological classification of Stokes complexes 
for polynomials of degree 3 is given. 

In the rest of the paper we apply these results to problems with PT- 
symmetry. In Section 4, we consider the PT-symmetric cubic family 
( II. 4p with real a and A. We prove that the intersection of the spectral 
locus with the real (a, A)-plane consists of disjoint non-singular analytic 
curves r„, n > 0, the fact previously known from numerical compu- 
tation [T^ H3| [3Tj . Moreover, we prove that the eigenfunctions corre- 
sponding to (a. A) G r„ have exactly 2n zeros outside the imaginary 
line. (They have infinitely many zeros on the imaginary line). Further- 
more, using the result of Shin on reality of eigenvalues, we study the 
shape and relative location of these curves r„ in the (a, A)-plane and 
show that a — )■ -|-oo on both ends of r„, and that for a > 0, F^ con- 
sists of graphs of two functions, that lie below the graphs of functions 
constituting F„+i. 

This gives PT-analog of the familiar fact for Hermitian boundary 
value problems that "n-th eigenfunction has n real zeros" ; in our case 
we count zeros belonging to a certain well-defined set in the complex 
plane. This result proves rigorously what can be seen in numerical 
computations of zeros of eigenfunctions by Bender, Boettcher and Sav- 
age [g. 

The result of Section 4 also gives a contribution to a problem raised 
by Hellerstein and Rossi |2]: describe the differential equations 

(1.6) y" + Py = 

with polynomial coefficient P which have a solution whose all zeros are 
real. For polynomials of degree 3, all such equations are parametrized 
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by our curve Fq, and equations having solutions with exactly 2n non- 
real zeros are parametrized by F^. 

The arguments in Section 4 use our parametrization of the spectral 
loci from [211 CH] combined with the singular perturbation results of 
Sections 2 and 3. These perturbation results allow us to degenerate the 
cubic potential to a quadratic one (harmonic oscillator) and to make 
topological conclusions based on the ordinary Sturm-Liouville theory. 

Next we apply similar methods to two families of PT-symmetric 
quartics 

(1.7) - w" + {z^ + az^ + icz)w = \w , m;(±oo) = 0. 
and 

(1.8) w" + {z^ + 2az^ + 2imz + X)w = 0, 

(1.9) w{re^) ^ 0, as r -^ oo, ^ e {-7r/6, -57r/6}, 

where m > 1 is an integer. The first family was considered in [3] and 
[121 [I3|- We prove that the spectral locus in the real (a, c, A)-space 
R^ consists of infinitely many smooth analytic surfaces Sn, n > 0, 
each homeomorphic to a punctured disc, and that an eigenfunction 
corresponding to a point (a, c, A) G Sn has exactly 2n zeros which do 
not lie on the imaginary axis. We study the shape and position of these 
surfaces by degenerating the quartic potential to the previously studied 
PT-symmetric cubic oscillator. 

The second quartic family (1 1.8m. 91) was introduced by Bender and 
Boettcher [1]. It is quasi-exactly solvable (QES) in the sense that 
for every integer m > 1 in the potential, there are m "elementary" 
eigenfunctions, each having m — 1 zeros. The part Zm of the spectral 
locus corresponding to these elementary eigenfunctions is a smooth 
connected curve in C^ [191 120] • In the end of Section 4 we study the 
intersection of this curve with the real (a, A)-plane. Similarly to the case 
of the PT-symmetric cubic, this intersection consists of smooth analytic 
curves Fj^^, n = 0, . . . , [m/2], and for (a. A) G T^^ the eigenfunction 
has exactly 2n zeros that do not lie on the imaginary axis. For n < m/2 
the part of T^ „ over some ray a > am consists of disjoint graphs of two 
functions, and we have the following ordering: (a. A) G FJ„„, (a. A') G 
^mn+i ^^^ 0, > ttm imply that A' > A. Moreover, the QES spectrum 
for a > am consists of the m smallest real eigenvalues. 

The results of Section 4 permit us to answer the question of Heller- 
stein and Rossi stated above for polynomial potentials of degree 4: 
All equations (II. 6p that possess a solution with 2n non-real zeros are 
parametrized by our curves FJ„„ if the total number of zeros is m — 1, 
and by our surfaces Sn if the total number of zeros is infinite. 
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Notation and conventions. 

1. What we call Stokes lines is called by some authors "anti-Stokes 
lines" and vice versa. We follow terminology of Evgrafov and Fedoryuk 
[231121]. 

2. We prefer to replace z by iz in PT-symmetric problems. Then 
potentials become real, and the difference between PT-symmetric and 
self-adjoint problems is that in PT-symmetric problems the complex 
conjugation interchanges the two boundary conditions, while in self- 
adjoint problems both boundary conditions remain fixed by the sym- 
metry. The main advantage for us in this change of the variable is 
linguistic: we frequently refer to "non-real" zeros. The expression 
"non-real" excludes 0, while the expression "non- imaginary" does not. 

We thank Kwang-Cheul Shin for his useful remarks and for sending 
us the text of his lecture [39] and Per Alexandersson for the plots of 
Stokes complexes he made for this paper. 

2. Perturbation of eigenvalues and eigenfunctions 

We begin with recalling some facts about boundary value problem 
(11.11) with potential (11. 2p . The separation rays are defined by 

Re ( / ^/c^ dn = 0, that is z'^+2 < 0. 

These rays divide the plane into d + 2 open sectors Sj which we call 
Stokes sectors. We enumerate them by residues modulo d + 2 counter- 
clockwise. 

A solution w of the differential equation (II. ip is called suhdominant 
in Sj if w{rz) — )■ as r — )■ +oo, for all z G 5*^. For every j, the space of 
solutions of the equation in (II. ip which are subdominant in Sj is one- 
dimensional. If Sj and Sk are adjacent, that is j = /c±l mod {d+2) then 
the corresponding subdominant solutions are linearly independent. 

Let Sj and Sk be two non-adjacent Stokes sectors. We consider the 
boundary conditions 

(2.1) w is subdominant in Sj and Sk- 

Such boundary value problem has an infinite set of eigenvalues tend- 
ing to infinity. All eigenspaces are one-dimensional. These facts were 
proved by Sibuya [ID] whose main tool were special solutions normal- 
ized on one ray, which we call Sibuya solutions. Precise definition is 
given below. Our first goal is to prove continuous dependence of Sibuya 
solutions on parameters. 
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We consider a family of polynomial potentials with parameters (t, a): 

m 

(2.2) Q{z,t,ai) = tz'^ + y^ajZ^ , m < d, a = (ao, . . . , 0^)- 

Let K C C"^+'^ be a compact set which has a fundamental system of 
open simply connected neighborhoods, and such that a^ 7^ for a G K. 
This compact K will be fixed in all our arguments, so our notation does 
not reflect dependence of the quantities introduced below on K. Let 
L = {te*^ : t > to} be a ray in C. 

Suppose that for some 6 > and for all (t, a) G [0, 1] x K the 
following conditions hold: 

a) There exists R > such that | arg2 — 6\ > S for all zeros z of 
Q{z,t,a.) such that \z\ > R. 

b) Whenever L intersects a vertical trajectory of Q{z,t,Si)dz'^ at a 
point z, \z\ > R, the smallest angle between this trajectory and L at 
the intersection point is at least 6. 

c) 6 is not a separating direction of Q(z, t, a.)dz^. 

One can easily show that b) implies c). Condition c) simply means 
that L is neither a Stokes line for tz'^sz'^, nor a Stokes line for a^z'^dz'^. 

Condition a) implies that there is a branch Ql of Q^ analytic on 
[0, 1] X /T X L. Let Ql = (Q^ )^ be the corresponding branch of Q^^^. 
We choose the original branch Q/ in such a way that 

ReQl {z)dz ^ +00, z — 7- 00, z E L. 

This is possible in view of condition c). 

Let us say that y = yL{z, t, a) is a Sibuya solution of 

-y" + Q{z,t,Si)y = 

if 

y{z)r^Ql'^\z)exp(-jy]!\OdcY z^oo, zEL. 

Here zq = toe^^. Notice that a change of to, results is multiplying the 
Sibuya solution by a factor that depends only on t and a but not on z. 

Theorem 2.1. Under the conditions a), b), c) above, there exists a 
unique Sibuya solution. It is an analytic function of {z, t, a) in a neigh- 
borhood ofCx(0,l]xK, continuous on Kix[0,l]xK for every compact 
Ki C C 
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Proof. Let 



J zo 



where the integral is taken along L This is an analytic function which 
maps L onto a curve 7. Condition b) implies that 7 is a graph of a 
function (intersects every vertical line at most once), the slope of 7 
is bounded, and maps bijectively some neighborhood of L onto a 
neighborhood of 7. 

Let ip be the inverse function to 0. 

Setting u = {Q/ o ijj^y o ■?/;, we obtain the differential equation 

where the primes stand for differentiation with respect to (, and 



/2 



5Q'' 1Q"\ 

This is equivalent to the integral equation 

(e^-^' - 6"^'+") g{r])u{r])dr], 

where the path of integration is the part of 7 from C to 00. The integral 
equation is solved by successive approximation. We set u{C) = f (C)e^, 
and obtain 



1 r°° 

(2.4) v{C) = 1 + - / (e^^'^"^) - 1) g{vHv)dv 



2-/C 



which we abbreviate as f = 1 + F{v). Setting f = and Vn+i 
1 + F{vn), we obtain 



\\Vn+l - VnWoo < \\Vn " ^'n-l|| / \9it)\dt. 

Here we used the fact that 5R(C — 77) < on the curve of integration 
because 7 is a graph of a function. Now, if C > is large enough, we 
have 

/oo 
\giv)\dv<l/2 

for all values of parameters t G [0, to], a G K. We state this as a lemma: 

Lemma 2.2. There exists h E L such that for the piece Lj, of L from 
b to 00 and for all (t, a) G [0, 1] x K we have 




2 



1" 



1(9^ _^_Q_]Q-y^dz 



< 1/2. 
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The integral in this lemma equals to the integral in (12 .Sp by the 
change of the variable z = V'(C)) VQdz = d(. 

Proof. Let ^i, . . . ,2;^ be all zeros of Qt listed with multiplicity, in 
an order of non-decreasing moduli. Suppose that zi, . . . , zm are in the 
disc \z\ < R while the rest are outside. Here R is the number from 
condition a). We have 

^ = 1:^+ E ^ = a,iz) + a,iz), 
and 

First we estimate 

(2.6) l^il < -, -, \a[{z)\ < 



\z\-R' ' ^^ ^' - (|z|-i?)2 
To estimate (T2 we first use condition a) to conclude that 

(2.7) \z - Zk\ > Co\z\, z e L, k>M 
where Cq depends only on 6. Then 

(2.8) \a2{z)\<C^/\z\, \a',\ < C2/\z\^ 
Applying these inequalities, we obtain that 



5 Qf 1 Q': 



- z^ 



16 g? 4 Q^ 

on L, where C depends only on 6. 
Now we write 

M d d 

\Qt{z)\ = tl[\z - z,\ H \z-z^\>tC',{\z\-R\)^' n l^il- 

j=l j=M+l j=M+l 

Here we used inequality (12. 7p with interchanged z and Zk- It is easy to 
see by Vieta's theorem that 

d 

n \zj\>CAt-\ 

j=M+l 

where C4 depends only on K. This shows that |Qt(-2)| > Cgl^l*'^. 
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Now we use the fact that L;, is a part of a ray from the origin, so 
\dz\ = d\z\ on Lh Putting all this together we conclude that our integral 
is majorized by the integral 



-*-^/2rf|; 



which proves the lemma. 



So the series ^fn is convergent uniformly in Re^ > A for some 
A > and this convergence is uniform with respect to t, a. Then an 
application of the theorem on uniqueness and continuous dependence 
of initial conditions for linear differential equations shows that this 
convergence is uniform also on compacts in the right half-plane of the 
C-plane. D 

Let Zt be a family of discrete subsets of the complex plane, depending 
on a parameter t. We say that Zt depends continuously on t if there 
exists a family of entire functions ft ^ such that Zt is the set of 
zeros of ft, and ft depends continuously on t. Here the topology on the 
set of entire functions is the usual topology of uniform convergence on 
compact subset of the complex plane. 

Consider the eigenvalue problem 

(2.9) -y" + Qy = Xy, y{z) ^ 0, z ^ oo z G Li U La, 

where Q is a polynomial in z of the form (12.21) . and Li,L2 are two 
rays from the origin in the complex plane. We say that a ray L is 
admissible if conditions a), b) and c) in the beginning of this section 
are satisfied. The notion of admissibility depends on the parameter 
region K participating in conditions a), b) and c). 

Theorem 2.3. // both rays Li and L2 are admissible then the spectrum 
of problem Ii2.9\) is continuous for (t, a) G [0, 1] x K. 

Proof. Let yi and 1/2 be the Sibuya solutions corresponding to the rays 
Li and L2. Then their Wronski determinant W = y[y2 — yiy'21 where 
the primes indicate differentiation with respect to 2, is the spectral 
determinant of the problem (12. 9p . and W depends continuously on 
(a. A) in view of Theorem 12.11 D 

The limit problem (12. 9 p for t = may have no eigenvalues, this is the 
case when the rays Li and L2 belong to adjacent sectors of Q{z, 0, a). 
In this case. Theorem 12.31 says that the eigenvalues escape to infinity 
as t — )■ 0. 
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Theorem 2.4. Consider the problem ^2. Pj) . and suppose that the rays 
Li and L2 are admissible, X{t, a) is an eigenvalue that depends con- 
tinuously on (t, a) and has a finite limit A(0,a) as t — )■ 0. Then there 
exists an eigenfunction y{z,t,a.) of this problem, corresponding to this 
eigenvalue, that depends continuously on (t,a). 



Proof. Let yi{z,t,a.) be the Sibuya solution of ( 12.9P with A = A(t,a), 
corresponding to the ray Li. It depends continuously on (t, a) by The- 
orem [231 Let 1/2(2, t, a) be the Sibuya solution of the same equation 
corresponding to L2. The assumption that X{z,t,a.)) is an eigenvalue 
implies that yi and ?/2 are proportional. This implies that yi tends to 
zero as 2; — )■ 00 on L2, so yi satisfies both boundary conditions. Thus 
yt is an eigenfunction that depends continuously on t. D 



3. Admissible rays 

In this section we give a criterion for a ray to be admissible (see the 
definition before Theorem 12. 3p . We reduce the problem to the case of 
a binomial Q{z) = tz'^ + cz'^, t ^ 0, c 7^ 0, < m< d. 

We begin with recalling terminology. Let Q{z) be a polynomial, 
2 e C. A vertical line of Q{z) dz"^ is an integral curve of the direction 
field Q{z) dz"^ < 0. A Stokes line of Q is a vertical line with one or both 
ends in the set of turning points {z : Q{z) = 0}. The Stokes complex 
of Q is the union of the Stokes lines and turning points. Examples of 
Stokes complexes are shown in Figs. [T]l3l 





4 I ,-^3 



(a) z + iz 



(b) z^ + e'^'I'^z^ 



Figure 1. Stokes complexes of z^ + iz^ and z'^ + e^'^/'^z^. 



A horizontal line of Q is a vertical line of —Q. Vertical and horizontal 
lines intersect orthogonally. An anti-Stokes line of Q is a Stokes line of 
-Q. 
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(a) z"* + z^ (b) z^ - z^ 

Figure 2. Stokes complexes of z^ + z^ and z^ - z^ . 




(a) z3 + ^2 (b) 




IZ 



Figure 3. Stokes complexes of z^ + z^ and z^ + iz} . 

Every Stokes line has one end at a turning point and the other end 
either at a different turning point or at infinity. If Q has a zero at Zq 
of multiplicity m then there are ?ti + 2 Stokes lines with the endpoint 
at Zi^\ they partition a neighborhood of Zq into sectors of equal opening 
271 /{m + 2). The m + 2 anti-Stokes lines having one end at Zq bisect 
these sectors. 

Let L(a) = {2 G C \ {0} : arg2; = a} and, for < /3 — a < 27r, 
^(a, /3) = {2 G C \ {0} : a < arg 2 < /?}. For i? > 0, let L>(i?) = {z e 
C : \z\ < R}. For a ray L or a sector S, define Lr = L\ D{R), Sr = 
S \ D{R). For a set 5 C C, :S is its closure in C and 9^ = ;S \ S. 

Definition 3.1. Let Q = Pd+Pm be a binomial, where Pd{z) = tz'^ and 
Pm{z) = cz"^ are two non-zero monomials, < m < d. Let S{Q) be 
the partition of C\{0} into open sectors and rays defined by the Stokes 
lines of the two monomials P^ and Pm- Let TZi^Q) be the refinement of 
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S{Q) defined by the rays from tlie origin tlirougli tlie non-zero turning 
points of Q. 

A ray L{6) is called good for Q if it is not one of the rays of TZiQ) and 
is not tangent to any vertical line of Q. The last condition is equivalent 
to L{9) n Z = where Z = {z : z^Q{z) < 0}. A sector S of n{Q) is 
good for Q if each ray L{6) C S* is good. This is equivalent to S'flZ = 0. 

Theorem 3.2. Let Q = Pd + Pm be as in Definition {3Ai Then any 
sector ofTZ{Q) containing an anti-Stokes line of Pd is good for Q, and 
any good ray belongs to one of such sectors. 

Proof. Let R_|_ and R_ be the positive and negative real rays (not 
including 0). Definition 13.11 implies that a ray L = L{6) is good if and 
only if the cone Cl = {az^Pd{z)+(5z^Pm{z), z e L, a e R+, /3 G R+} 
does not contain R_. 

An anti-Stokes line of Pd is a good ray unless it is also a Stokes line 
of Pm, because z^Pd{z) is real positive on it, hence z'^Pm{z) must be 
real negative to make the sum real negative. 

Suppose that L is an anti-Stokes line of Pd, and that z'^Pm{z) is either 
real positive or belongs to the upper half-plane for z & L. Then either 
Cl = R+ or Cl belongs to the upper half-plane. When L is rotated 
counterclockwise, the arguments of the two monomials z'^Pd and z^Pm 
restricted to L are increasing. Hence Cl remains in the upper half- 
plane until at least one of the monomials becomes real negative on L, 
i.e., until L becomes a Stokes line of either Pd or P^., or both. When L 
is rotated clockwise, the arguments of the two monomials restricted to 
L are decreasing. Since the argument of Pd decreases faster than the 
argument of Pm, the cone Cl does not contain negative real numbers 
until either z^Pd on L becomes real negative or argP^ — argP^ passes 
TT, i.e., until L either becomes a Stokes line of Pd or passes a non-zero 
turning point of Q. 

The case when Pm{z) belongs to the lower half-plane on an anti- 
Stokes line of Pd is done similarly. 

Conversely, let L be a ray which is not one of the rays of 'R-{Q) and 
such that Cl does not contain negative real numbers. Then either L 
itself is an anti-Stokes line of Pd, or it can be rotated to the closest 
anti-Stokes line of Pd preserving this property. 

For example, if z^Pd{z) belongs to the upper half-plane for z E L, 
then aig^z"^ Pd{z)) — vr < aigi^z"^ Pm{z)) < vr for 2 G L. Otherwise, either 
L would be a Stokes line of Pm (if arg(2;^Pm(2)) = vr), or it would 
contain a turning point of Q (if a.Tg{z'^Pd{z)) — 11 = aigi^z"^ Pm{z))) , 
or Cl would contain R_. When L is rotated clockwise, since argP^ 
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decreases faster than argP^, Pd would become real positive on L before 
either argP,„ — argP^^ becomes it or z'^P^ becomes real negative. 

The case when z'^P^iz) belongs to the lower half-plane on L is done 
similarly, rotating L counterclockwise. D 

Corollary 3.3. Let S = S{a,P) be a Stokes sector of Pa such that the 
anti- Stokes line L G S of Pd does not contain a turning point of Q. 
Then S contains a good subsector. 

Proof. Since L does not contain a turning point of Q, it cannot be a 
Stokes line of Pm- Hence L belongs to a sector of IZi^Q), which is good 
by Theorem 13. 2[ D 

Theorem 3.4. Let Q = Pd + Pm be as in Definition \3.1[ Let H = 
{{x, y) G R^ : |x| < vr, \y\ < tt , \y — x\ < it} be a hexagon in R^. Then 
R4. is a good ray for Q if and only i/(argt, argc) G H. Here the values 
o/argt and argc are taken in (— 7r,7r]. 

A ray L{6) is a good ray for Q if and only if (argt, argc) belongs to 
H translated by (27r/c — {d + 2)9, 2nl — (m + 2)6') for some integers k 
and I. 

Proof. For t > and c > 0, R+ is an anti-Stokes line of both Pd and 
Pm, hence it is a good ray. It is not a Stokes line of either Pd or Pm 
when I argt| < vr and | argc| < vr. It does not contain a turning point 
if I argt — argc| < vr. 

For I argt| < ir, the anti-Stokes line L of Pd closest to R+ has the 
argument —a.Tgt/{d + 2). It is a Stokes line of Pm if argc — (m + 
2) argt/((i + 2) = 'iT{2k + 1) for an integer k. Since < m < d, the 
lines y — {m + 2)x/{d + 2) = 7r(2/c -|- 1) do not intersect H. Hence L 
does not contain a turning point of Q when (argt, argc) G H. This 
implies that the closure of the sector 5* bounded by R4. and L does not 
contain Stokes lines of either Pd or Pm, and does not contain non-zero 
turning points of Q, for all (t, c) such that (argt, argc) G H. Due to 
Theorem 13.21 R+ is a good ray for Q with these values of t and c. 

If one of the three inequalities defining H becomes an equality, R+ 
becomes either a Stokes line of one of the two monomials of Q, or 
contains a turning point of Q. 

When (argt, argc) crosses one of the two segments |argt| < vr, 
I arg c| < vr, I arg t — arg c| = vr of the boundary of H, a turning point of 
Q crosses R+ and remains inside the sector S for all values (t, c) such 
that I argt| < TT, |argc|<7r, | argc — (m + 2) argt/((i + 2)| < vr. Due 
to Theorem 13. 2[ R_|_ is not a good ray for Q with these values of t and 



c. 



This implies that R+ is not a good ray for Q when (argt, argc) ^ H. 
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The statement for L{6) is reduced to the statement for R+ by the 
change of variable z = ue^^ in the quadratic differential Q{z) dz^. D 

Example 3.5. Let us investigate when the two rays of the real axis 
are good for a binomial Q{z) = tz'^ + C2™. According to Theorem 13. 4[ 
R+ is a good ray when (argt, argc) G H (green hexagon in Fig. H^). 



71 



-71 



-71 



71 




{a) 



(b) 





(c) 



id) 



Figure 4. The values of (argt, argc) in Example 13.51 
where both R_|_ and R_ are good (green) and where R_ 
is not good (red), (a) d and m even; (b) d and m odd; 
(c) d even, m odd; (d) d odd, m even. 



For R_ (6* = vr in Theorem 13. 4p there are 4 cases, depending on the 
parity of d and m. 

a) If d and m are even, both R+ and R_ are good rays for any Q 
such that (argt,argc) G H. 

b) If d and m are odd, R_ is a good ray for Q when (arg t, arg c) 
belongs to the complement in {—it, tt]^ of the two triangles (red area in 
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Fig.gb) with the vertices (0, 0), (0, tt), (-vr, 0) and (0, 0), (tt, 0), (0, -tt), 
respectively. Both R+ and R_ are good rays for Q when (arg(t), arg(c)) 
belongs to the union of the two squares (0, 1)^ and (—1, 0)^ (green area 
in Fig. Hb). 

c) If d is even and m is odd, R_ is a good ray for Q when (arg t, arg c) 
belongs to the complement in (— vr, vr] of the two triangles (red area 
in Fig. St) with the vertices (0,0), (vr,0), (vr,vr) and (0,0), (-vr,0), 
(— vr, — vr), respectively. Both i?+ and i?_ are good rays for Q when 
(arg(t), arg(c)) belongs to the union of the two open parallelograms 
(green area in Fig. Ht) with the generators (tt, 0), (— vr, — vr) and (vr, vr), 
(— vr, 0), respectively. 

d) If d is odd and m is even, R_ is a good ray for Q when (arg t, arg c) 
belongs to the complement in (— vr, vr]^ of the two triangles (red area 
in Fig. Hi) with the vertices (0,0), (0,vr), (vr, vr) and (0,0), (0, -vr), 
(— vr, — vr), respectively. Both /?+ and i?_ are good rays for Q when 
(arg(t), arg(c)) in the union of the two open parallelograms (green area 
in Fig. mi) with the generators (0,vr), (— vr, — vr) and (0, — vr), (vr,vr), 
respectively. 

In particular, if t is on the positive imaginary axis, both i?_|_ and /?_ 
are good rays for Q when — vr/2 < arg(c) < vr in case (a), < arg(c) < vr 
in case (b), —vr/2 < arg(c) < or vr/2 < arg(c) < vr in case (c), 
—vr/2 < arg(c) < vr/2 in case (d). 

By definition, a good ray L is not tangent to the vertical lines of Q 
and is not a Stokes line of either P^ or P^- Since the angles between 
L and vertical lines of Q have non-zero limits at the origin and at 
infinity, there is a lower bound for these angles on L. This lower bound 
depends continuously on L, hence there is a common lower bound for 
these angles for all rays in a proper subsector T of a good sector S 
(such that T \ {0} C S). 

The good sectors in Theorem 13.21 depend continuously on the argu- 
ments of the coefficients t and c of monomials Pd and Pm, except when 
a good sector degenerates to a ray that is a Stokes line of Pm and an 
anti-Stokes line of P^. 

The lower bounds for the angles between a good ray L and vertical 
lines, and for the values of P, depend continuously on the arguments 
of t and c, except when a good sector containing L degenerates. 

Now we show that a good ray for a monomial tz'^ + cz"^ is admissible 
for the potential (12. 2p . with 

K = {{ao,...,am-i,c) : sup \aj\<M}, 

0<j<m-l 

for every positive M. 
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Lemma 3.6. Consider the polynomials Q{z) = tz'^ + cz^ + q{z) , where 
^ > 0, \q{z)\ < M\zY^~^ and \c\ = 1. Let S be a sector whose closure 
does not contain turning points of tz'^ + cz"^ . 

For every e > there exists R > depending on e, M, S, such that 
for every t > 0: 

(i) The set S r\ {z : \z\ > R\ does not contain turning points of Q, 
and 

(a) Ifv{z) andv\z) are the vertical directions ofQ[z)dz^ and {tz'^ + 
cz"^)dz^ , respectively, then \ 8JLgv{z)—a,Tgv'{z)\ < e for z G S, \z\ > R. 

Proof We have \tz'^ + cz'^\ > cj^p for z e S,so \q{z)/{tz'^ + cz"")] < 
'^, and (i), (ii) hold when R is large enough. D 



c ^\z 



4. PT-SYMMETRIC POTENTIALS AND LINEAR DIFFERENTIAL 
EQUATIONS HAVING SOLUTIONS WITH PRESCRIBED NUMBER OF 

NON-REAL ZEROS 

Hellerstein and Rossi asked the following question [9l Problem 2.71]. 
Let 

(4.1) w" + Pw = 

be a linear differential equation with polynomial coefficient P. Char- 
acterize all polynomials P such that the differential equation admits a 
solution with infinitely many zeros, all of them real. 

This problem was investigated in [IS [291 [28l [33l [381 [22] . Recently 
K. Shin [39] announced a description of polynomials P of degree 3 or 4 
such that equation (14.11) has a solution with infinitely many zeros, all 
but finitely many of them real. It turns out that equations (14. ip with 
this property are equivalent to (11.41) or (11.71) of the Introduction by an 
affine change of the independent variable. 

Here we use the methods of [2T1 [TH] to parametrize polynomials P of 
degrees 3 and 4 such that equation (14.11) has a solution with prescribed 
number of non-real zeros. 

We begin with degree 3. 

Theorem 4.1. For each integer n > there exists a simple curve r„ 
in the plane R^ which is the image of a proper analytic embedding of 
the real line and which has the following properties. 
For every (a. A) G r„ the equation 

(4.2) - w" + {z^ -az + X)w = 

has a solution w with 2n non-real zeros. Real zeros belong to a ray 
(— cx),a;o) and there are infinitely many of them. This solution satisfies 
\imt^±^w{it) = 0. 
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The union U^o^n coincides with the real part of the spectral locus 

The projection (a, A) i— ?■ a, 

r„ n {(a, A) : a > 0} ^ {a : a > 0} 

is a 2-to-l covering map. The curves r„ are disjoint, and for a > 
and n > 0, if (a, A) G r.„ and (a, A') G r.„+i then A < A'. 




Figure 5. Curves T^ 
(Trinh, 2002). 



n 



0, . . . , 4 in the (a, A) plane 



Equation (14. 2p is equivalent to the PT-symmetric equation fll.4p in 
the Introduction by the change of the independent variable z i— )■ iz. 
Computer experiments strongly suggest that the projection (a, A) i— )• 
a is 2-to-l on the whole curve r„ except one critical point of this 
projection, and that the whole curve F^+i lies above r„. 

Fig.O taken from Trinh's thesis |13] (see also [S]), shows a computer 
generated picture of the curves F„. 

As a corollary from Theorem 14.11 we obtain that every eigenvalue 
A„(a), a > of (11. 4p . when analytically continued to the left along the 
a-axis, encounters a singularity for some a < 0. According to Theorem 
2 of [19] this singularity is an algebraic ramification point. 
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Proof. Consider the Stokes sectors of equation (14.21) . We enumerate 
them counter-clockwise as S'o, . . . 5*4 where 5*0 is bisected by the positive 
real axis. Consider the set G of all real meromorphic functions / whose 
Schwarzian derivatives are real polynomials of the form —2z^+a2z'^+aQ, 
and whose asymptotic values in the sectors S'o, . . . , 5*4 are oo, 0, b,b, 0, 
respectively, where fe = e*^, /3 G (0, vr). Such functions are described by 
certain cell decompositions of the plane [19]. By a cell decomposition 
we understand a representation of a space X as a union of disjoint 
cells. This union is locally finite, and the boundary of each cell consists 
of cells of smaller dimension. The 0-cells are points, vertices of the 
decomposition. The 1-cells are embedded open intervals, the edges, and 
the 2-cells are embedded open discs, faces of a decomposition. Two cell 
decompositions of a space X are called equivalent if they correspond 
to each other via an orientation-preserving homeomorphism of X. 

To describe functions of the set G, we begin with the cell decompo- 
sition $ of the Riemann sphere shown in Fig. Ei 



Figure 6. Cell decomposition $ of the image sphere 
(solid lines). 

It consists of one vertex at 2 and three edges which are simple disjoint 
loops around b,b and oo, so that the loop around oo is symmetric 
with respect to complex conjugation while the loops about b and b 
are interchanged by the complex conjugation. The point is outside 
the Fig. [61 The dotted line is not discussed here; it is needed for the 
future. Also for the future use, we assume that the loop around oo 
passes through the point —2, and that this loop is symmetric with 
respect to the reflection z H- — ^. 

So our cell decomposition has one vertex, three edges and four faces. 
The faces are labeled by points b, b, oo and which are inside the faces. 
(So three faces are bounded by single edge each, while one face (labeled 
with 0) is bounded by three edges). 

Suppose now that we have a local homeomorphism g : C ^ C such 
that the restriction 

(4.3) g:C\g~\A)^C\A, 



20 



ALEXANDRE EREMENKO AND ANDREI GABRIELOV 



where A = {6,6, oo,0}, is a covering map, and ^E' = g^^{^). Then 
we preimage \1/ = g~^{^) will be a cell decomposition of the plane 
C. Now suppose that a cell decomposition \E' of the plane is given 
in advance, and suppose that its local structure is the same as that 
of $. This means that the faces of "$ are labeled by elements of the 
set A, and that a neighborhood of each vertex of ^ can be mapped 
onto a neighborhood of the vertex of $ by an orientation-preserving 
homeomorphism, respecting the labels of the faces. Then there exists 
a local homeomorphism (7 : C — )■ C such that (14. 3p is a covering map. 
We use the following cell decomposition to construct g: 




Figure 7. Cell decomposition \E'„ for n = 2 (solid lines). 
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The five "ends" extend to infinity periodically. This cell decompo- 
sition depends on one integer parameter n > which is the number 
of 0-labeled faces between the neighboring "ramification points" . Only 
some face labels are shown but the reader can easily recover all other 
labels from the condition that in a neighborhood of each vertex \i/„ 
is similar to a neighborhood of the vertex of $. The dotted lines are 
not a part of our cell decomposition; they are preimages of the dotted 
line in Fig. 19, and are added for a future need. ^„ is symmetric with 
respect to the real line, this permits to make our local homeomorphism 
g symmetric, that is g{'z) = g{z). This construction defines the map g 
up to pre-composition with a symmetric homeomorphism : C — ?■ C 
of the domain oi g. A fundamental result of R. Nevanlinna ensures 
that this homeomorphism can be chosen in such a way that f = go4> 
is a meromorphic function which is real in the sense that f{z) = f{z). 
We refer to [19] for the discussion of this construction in our current 
context; in fact [I9] contains a simple alternative proof of Nevanlinna's 
theorem. Nevanlinna's original proof is explained in modern language 
in [TS]; the original paper of Nevanlinna is [SB] . 

The meromorphic function / is defined by the cell decomposition \i/„ 
and parameter b up to pre-composition with a real affine map cz + d. 
Furthermore, the Nevanlinna theory says that the Schwarzian deriva- 
tive of / is a polynomial of degree exactly 3 (the number of unbounded 
faces of "^n minus 2). We pre-compose / with a real afiine map to 
normalize this polynomial to have leading coefficient —2 and zero co- 
efficient at z^. Thus 

(4.4) ^"2(7) =-2(^'-«^ + ^)- 

As / is real, a and A are also real. Now / is uniquely defined by 
the properties that it satisfies a differential equation (14. 4p . has asymp- 
totic values 00, 0, b,b, in the sectors 5*0, ... , 84, respectively, and that 
/~^($) equivalent to \l/„ (Fig. [7]for n = 2) by an orientation-preserving 
homeomorphism of the plane commuting with the reflection 2; 1— )■ z. 
The statement on asymptotic values implies that 

f{it) ^0, t -> ±00. 

Furthermore, / depends analytically on b, when b is in the upper half- 
plane, and thus we obtain a real analytic map b H- (a. A). This map is 
evidently invariant with respect to transformations b \-^ tb, t G R\{0}, 
this is because the Schwarzian derivative in the right hand side of (14. 4p 
does not change when / is replaced by tf. 
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Thus for every n, we have a one-parametric family G„ C G of mero- 
morphic functions, parametrized by /3 G (0, vr), 6 = e^^ . Taking the 
Schwarzian derivative we obtain a map F„ : (0,7r) — > R^, /3 i— ;■ (a, A). 
This map is known to be a proper real analytic immersion [2]. It is 
easy to see that it is injective: two solutions of the same Schwarz equa- 
tion may differ only by post- composition with a fractional-linear map, 
and this fractional-linear map must be identity by our normalization 
of asymptotic values. 

For the same reasons the images of F^ are disjoint: for different n, 
our functions have (topologically) different line complexes. The images 
of Fn are our curves r„. 

Now we prove that the union of F^ coincides with the real part of 
the spectral locus of (11. 7p . 

Our functions f & G can be written in the form / = w/wi where w 
and Wi are two linearly independent solutions of equation (14. 2 p with 
some real a and A. We can choose w and Wi to be real entire functions. 
Condition that f{it) — )■ as t — )■ ±oo implies that w(2t) — )■ for 
t — )■ ±oo so w is an eigenfunction of the spectral problem 

(4.5) -w" + {z'^ -az + \)w = Q, w(±20o) = 0, 

which is equivalent to (11.41) by the change of the independent variable 
z i-> iz. 

So our curves F„ belong to the real part of the spectral locus of (14.51) 
or ([HD. 

Now, let A be a real eigenvalue of the problem (14. 5p . w a corre- 
sponding eigenfunction. Choose a point xq on the real line such that 
w{xq) 7^ and normalize w so that w{xq) = 1. Then w*{z) = w(z) 
is an eigenfunction with the same eigenvalue, so w* = cw for some 
constant c 7^ 0. Substituting xq gives that c = 1. So w is real. 

Let wi be a solution of the same equation as w but satisfying w{x) -^ 
0, a; — )■ -|-oo. We normalize Wi so that Wi is real in the same way as we 
normalized w. Then / = w/wi is a real meromorphic function whose 
Schwarzian derivative is a cubic polynomial with top coefficient —2, 
and the asymptotic values in Sj are oo,0,b,b,0. We can change the 
normalization of Wi multiplying it by any real non-zero constant. In 
this way we achieve that b = e*^ for some /3 G (0,7r). So / belongs to 
the class G. 

Lemma 4.2. G = U^^oGn- 

Proof. Let f E G. Consider the cell decomposition X = /~^($). We 
have to prove that X = \E'„ for some n > 0. To do this, we follow 
[19]. We first remove all loops from X, and then replace each multiple 
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edge by a single edge, and denote the resulting cell decomposition by 
Y. Notice that the cyclic order (00,6,6) in Fig. [6] is consistent with 
the cyclic order (00, 0,6, 6,0) of the Stokes sectors in the z-plane. By 
[T9| Proposition 6] , this implies that the 1-skeleton of F is a tree. This 
infinite tree is properly embedded in the plane, has 5 faces, is symmetric 
with respect to the real line, and has two faces labeled with which are 
interchanged by the symmetry. Moreover, the faces of Y are in one-to- 
one correspondence with the Stokes sectors, and the face corresponding 
to 5*0 is bisected by the positive ray. One can easily classify all trees 
with these properties. They depend of one integer parameter n > 
which is the distance between the ramification point in the upper half- 
plane and the ramification point on the real axis. Now we refer to [191 
Proposition 7] that the tree Y uniquely defines the cell decomposition 
X. This shows that X = \Ef„ for some n > 0. 

Meromorphic function / is defined by the cell decomposition X and 
the parameter 6 up to an afiine change of the independent variable. 
Normalizing it as in fl4.4p gives / G (?„. D 

We conclude that the union of our curves r„ in the right half-plane 
a > coincides with the real part of the spectral locus of f l4.5p . 

Now we study the shape of the curves F^. The boundary value 
problem fl4.5p was considered by Shin [37] , Delabaere and Trinh ^\ [T3] . 
The spectrum of this problem is discrete, simple and infinite. It is 
known [3^ that for a > all eigenvalues of this problem are real and 
positive. It follows from this result that there are real analytic curves 
A = 7fc(a), k = 0,1,2, ... , such that for each k, 7fc(a) is an eigenvalue 
of the problem (14. 5p . and 7^ < 7^+1, A; = 0, 1, 2, . . .. So the part of the 
real spectral locus in {(a. A) : a > 0} is the union of the graphs of 7^. 

Next we prove that the intersection of F„ with the half-plane a > 
consists of 72n and 72n+i- For this purpose we study what happens to 
eigenvalues and eigenf unctions of the problem (14. 5 p as a ^^ +00. 

A different approach to the asymptotics as a — )• 00 is used in [27] . 
We could use their Corollary 2.16 here instead of referring to Sections 
2,3. 

We substitute cz + d in (14. 5 p and put y{z) = w{cz + d), where 

rf=(a/3)i/2>o, c=(3ci)"^/^>0. 

The result is 

(4.6) -y" + {c^z^ + z^ + fi)y = 0, 

where /i = c^lX+d^—ad). Choosing the positive and negative imaginary 
rays as our normalization rays Li and L2, we see that the normalization 
rays are admissible in the sense of Theorem 12.31 The Stokes complex 
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of the binomial potential corresponding to (14. 6p is shown in Fig. |3]^a) . 
According to Theorem l2.3l the spectrum of the problem (14.61) converges 
to the spectrum of the limit problem 

(4.7) - y" + z^y = -jjy, y{±ioo) = 0. 

This limit problem is equivalent to the self-adjoint problem 

—u" + z'^u = fiu, n(±oo) = 

by the change of the variable u{z) = y{iz). Convergence of the spec- 
trum implies convergence of eigenfunctions uniform on compact subsets 
of the plane by Theorem 12.41 As a varies from to oo, we can choose 
an eigenvalue A (a) which varies continuously, and the corresponding 
eigenfunction that varies continuously, and tends to an eigenfunction 
of (14. 6p . In the process of continuous change the number of non-real 
zeros of the eigenfunction cannot change because eigenfunctions cannot 
have multiple zeros. The conclusion of the theorem will now follow from 
the known properties of zeros of eigenfunctions of Hermitian boundary 
value problems, once we establish the following 

Lemma 4.3. As t = C' -^ Q in ^4-0^ the non-real zeros of an eigen- 
function cannot escape to infinity. 

Notice that the real zeros of the eigenfunction do escape to infinity, 
as the limit eigenfunction has at most one real zero. 

Proof. Let Wt be the eigenfunction constructed in Theorem 12.41 which 
depends continuously on t. Let w^ be the Sibuya solution corresponding 
to the positive ray. Then ft = Wt/wl is a real meromorphic solution 
of the Schwarz equation and has asymptotic values oo, 0, ht, h, in the 
sectors Sj. As ft — )■ /o, and the Schwarzian of /o is of degree 2, we 
conclude that b = e^'^ converges to the real axis, and the Riemann 
surface of /^"^ must converge in the sense of Caratheodory [101 SS], to 
a Riemann surface with 4 logarithmic branch points which can lie only 
over 0, oo,6o) where bo E {±1}. To construct the cell decomposition 
corresponding to this limit Riemann surface, we consider two cases. 

Case 1. 6o = 1- To describe the limit function, we must replace in 
the original cell decomposition Fig. |6] two loops corresponding to 6, b 
with a single loop around both of these points. This loop is shown by 
the dotted line in Fig. |6] and its preimage is shown by the dotted lines 
in Fig. [71 The original loops that separate b- and b- labeled faces from 
the face labeled must be removed. Performing this operation on the 
cell decomposition \l/„ we see that the 1-skeleton breaks into infinitely 
many pieces. But there is only one piece that has four unbounded faces 



SINGULAR PERTURBATION, PT-SYMMETRIC FAMILIES 



25 



and thus can correspond to a meromorphic function whose Schwarzian 
derivative is a polynomial of degree 2. This limit decomposition is 
shown in Fig. |8l 



■y.. 



"■I. 




h=h 



% 



A" 








Figure 8. Limit cell decomposition with n = 2 (solid 
and dotted lines). 

This time both solid and dotted lines represent the edges of this 
decomposition. We see that the number of non-real zeros in the limit 
is the same as it was before the limit. 

Case 2hQ = —\. To analyze this case, we replace the cell decompo- 
sition on Fig. |6] by the one in Fig. [9] 

Now we want to find the preimage of the cell decomposition in Fig. |9] 
under the same function /. There are several ways to find this preim- 
age. Let us choose for convenience b = i, and express the loops in 
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Figure 9. Another cell decomposition of the sphere 
(sohd lines). 



Fig. ini in terms of the loops in Fig. |6l as elements of the fundamental 
group of C\{0, i,—i}. We denote the loops around b,b in Fig. |6] by 
'^hilhi and let a, (3 be the upper and lower halves of the loop around 
oo, so that 7oo = a/3 (a followed by (3). Let 7^ and 7^ be the loops in 
Fig. ini Then we have 700 = 7^ = Pa (a followed by /3), 7^ = P'jbP'^, 



7 



a 



^Tba. See Fig. M 




Figure 10. Fig. [6] and Fig. [9] together. 



These relations permit us to draw the preimages of the loops 7^, 7^, 7^ 
in the z-plane. The resulting picture is shown in Fig. [HI 

When 6 — )■ — 1, we have a degeneration as before. The corresponding 
cell decomposition is obtained by replacing preimages of the loops 7^ 
and 7^ by the preimages of the dotted line. The resulting cell decom- 
position is shown in Fig. [121 

We see that the limit function still has 2n non-real zeros, and one 
real zero. This completes the proof of the lemma, and of Theorem 
EB □ 
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Figure 1 1 . Preimage of the cell decomposition Fig. [9] 
(solid lines). 



This proof shows that /3 — )■ corresponds to the lower branch of r„ 
while /9 — )■ TT corresponds to the upper branch. 

D 
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Figure 12. The limit cell decomposition of Fig. [TOlas 
6 — 7- — 1 (solid and dotted lines). 



Theorem 4.4. Let P he a polynomial of degree 3 such that equation 
( 1^. j| j has a solution with 2n non-real zeros. Then ( [^. j| j can be trans- 
formed to an equation of Theorem 4-1 with (a, A) G r,„ by a real affine 
change of the independent variable. 



Proof. By the results of Gundersen [2S1 ESj , all solutions have infinitely 
many zeros, and the coefficients of P are real. By a real affine change 
of the variable we achieve that P{z) = —z^ + az — \. As almost all 
zeros are real, our solution must tend to zero in both directions of the 
imaginary axis. So A is an eigenvalue of the problem (14. 5p . Let w 
be a real eigenfunction and Wi a real solution of our equation that is 
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linearly independent of w. Then the ratio / = w/wi is a meromorphic 
function which is a local homeomorphism, and has asymptotic values 
oo, 0, 6, 6, in So, . . . , S4, respectively. After a real afiine change of the 
independent variable, this function will belong to the class G defined 
in the proof of Theorem I4.1[ D 

Now we state analogous results for quartic oscillators. There are two 
different real two-parametric families in which solutions with finitely 
many non-real zeros can occur [32] ■ 

(4.8) - w" + {-z^ + az^ + CZ + X)w = 0, w(±ioo) = 0. 
studied in [31 |T3] , and 

(4.9) -w" + {z^-2az^ + 2mz + \)w = 0, w{te'^) -^ 0, ^ = ±7r/3 

studied in [1]. Here m > 1 is an integer. Problem 14.91 is quasi-exactly 
solvable, which means that there are m eigenfunctions of the form 
p{z)exp{z^/3 — az), where p is a polynomial of degree m — 1. The 

families (14. 8 p and (14. 9 p are equivalent to the PT symmetric families 
(ll.7p and (ll.8m.9p of the Introduction via the change of the independent 
variable z \-^ iz. 



Theorem 4.5. The real part of the spectral locus of ( [^.g[ ) consists 
of disjoint smooth connected analytic surfaces Sn, n > 0, properly 
embedded in R^. For (a, c. A) G Sn, the eigenfunction has 2n non-real 
zeros. Each of these surfaces is homeomorphic to a punctured disc. 
Projection 7r(a, c. A) = (a, c) has the following properties: It is a 2-to-l 
covering over some neighborhood of the a-axis, and for a > a^, the 
preimage of every line c = const is compact and homeomorphic to a 
circle. 

Proof. We follow the same pattern as in the proof of Theorem 14.11 
There are 6 Stokes sectors, Sq, . . . , S^, which we enumerate anticlock- 
wise, beginning from the sector in the first quadrant. 

If / = w/wi where w is a real eigenfunction and wi is a real linearly 
independent solution of the same equation, then / has asymptotic val- 
ues bo, 0, 61, bi, 0, bo in the sectors So, ... , S5. Here bo 7^ &i, and bo, 61 
must belong to C\R. 

If c = 0, we can choose w,wi with the additional symmetry with 
respect to the imaginary axis, which gives bo = —bi, so bo and 61 
belong to the same half-plane of C\R. The same situation persists for 
all real c because bo,bi depend continuously on c and never cross the 
real line. The real afiine group acts on / by post-composition; this 
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corresponds to the change of normahzation of w and Wi. So we can 
always choose the normahzation so that hi = i. 

Notice that after this normahzation condition c = corresponds to 
|6o — V2| = 1/2. See Remark [4.61 after the proof. 

Consider the ceU decomposition $ of the Riemann sphere (the range 
of /) shown in Fig. [131 It consists of one vertex at oo and four disjoint 
loops around ±2 and 6, h that are interchanged by the symmetry. 




Figure 13. Cell decomposition $ of the sphere (solid lines). 



Now consider the cell decomposition \l/„ of the plane (with labeled 
faces) shown in Fig. [Ml It is locally similar to $, and depends on one 
integer parameter n > which is the number of 0-labeled faces between 
the adjacent "ramification points". 

As in Theorem 14.11 Nevanlinna theory gives for each n > a family 
Gn of meromorphic functions / which have 2n non-real zeros and satisfy 
the Schwarz equation of the form 



(4,10) 



^-.M^ =2(. 



az 



cz 



A). 



f 2V/'. 

with real a, c, A. 

Classification result for symmetric trees with 6 faces in [12] ensures 
that all equations (14. ip having a solution with infinitely many real zeros 
and 2n non-real zeros are equivalent to equations which arise from our 
families G„. 

This also has an implication that there is "no monodromy" in our 
families Gn'- when h traverses a loop around i, we return with the 
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Figure 14. Cell decomposition \l/n for a quartic with 
n = 2 (solid lines). 



same function / we started with. Indeed, in the process of continuous 
deformation the number of non-real zeros cannot change, and there is 
only one suitable cell decomposition \l/„ for every n. 

Thus our family Gn is homeomorphic to a punctured disc. Taking the 
coefficients a, c, A of the Schwarzian defines an analytic embedding of 
Gn to R'^. This is our surface S'„. The surfaces are disjoint and properly 
embedded for the same reasons as in the proof of Theorem 14. 1[ 

To study the shape of these surfaces Sn in R'^, we first notice that 
for c = 0, the eigenvalue problem obtained from f l4.8p by rotation 
z ^-^ iz is Hermitian. It follows that the intersection with the plane 
Sn n {(a,c. A) : c = 0} consists of the disjoint graphs of two analytic 
functions defined for all real a, and that A.„(a, 0) < A„+i(a, 0). Another 
simple property of the surface Sn is that it is symmetric with respect 
to change c i— ;■ — c, which follows by changing z H- —z in the equation. 
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Now we study the asymptotic behavior of Sn for a — )■ +00. In the 
equation fl4.8p we set z = e{( — t), y{z) = w{e{( — t)), where t satisfies 

(4.11) a-6eV = 0, and AeH = 1, 
and obtain 

(4. 12) - y" + {-e^z^ + z^ + az + fj.)y = 0, 



where 
and 



a = Aeh^ + 2e^at + c£^ 



.6.4 I „^4,2 „^3. , ^2- 



/i = -e^r + aeH^ - ceH + e^X. 

Expressing t and a from equations (14.111) as functions of e and substi- 
tuting the result to the expression of a we obtain 

(4.13) a = (3/8)e-^°, c = e-^a - (l/4)e"^^ 

(4.14) A = -21 • 2-^20 + (a/4)5-8 + /i£-2. 

Consider the curves r„ from Theorem 14. 1[ It follows from their proper- 
ties stated in Theorem 14. H that for every n, there exists a„ = max{— a : 
(a, A) G r„, and < a„ < 00. 

Suppose that a < a„, and consider the curve in (a, c)-plane paramet- 
rized by (14.131) . Equation (I4.12p satisfies the conditions of Theorem 14. II 
of Section 2 (the Stokes complex corresponding to this equation is 



shown in Fig. 1(a), rotated by 90°), and the sectors containing the 
normalizing rays are stable. We conclude that the spectrum of (14.121) 
tends to the spectrum of the cubic 

(4.15) -y" + (z^ + az + fi)y = 0. 

The spectrum of the cubic with parameter a has at least one eigenvalue 
fi* which is real and such that the corresponding eigenfunction has 
2n non-real zeros. As fi* is an isolated point of this spectrum, and 
the spectrum of (I4.12p is symmetric with respect to the real axis, we 
conclude that there is an eigenvalue /i of (I4.12p which is real, and the 
corresponding eigenfunction has 2n non-real zeros. 

To ensure that the number of non-real zeros does not change in the 
limit, we make an argument similar to that in the proof of Theorem 
14.11 the degeneration of the cell decomposition \l/„ is shown in Fig. [T3 

We conclude that projection of our surface S contains a piece of the 
curve (I4.13P for e G (0,£:o)- 

Now suppose that a > «„. We claim that there are no points on Sn 
with a — )■ 00 and (a, c) on the curve (14.130 . Proving this by contra- 
diction, we suppose that there is a sequence (a^, Cj, Xj) G Sn such that 
{aj,Cj) belong to the curve (I4.13p . Then Theorem 12.31 implies that the 
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b=b 




Figure 15. Limit cell decomposition with n = 2 (solid 
and dotted lines). 



sequence /Xj related to the Xj by fl4.14p . has the property that /ij tends 
to a real eigenvalue fi* of the cubic oscillator fl4.15p . Then the corre- 
sponding eigenfunction tends to an eigenfunction of the cubic with 2n 
non-real zeros. This is a contradiction because a > a„, so our claim is 
proved. 

So the projection of S'„ on the plane (a, c) looks as a paraboloid 
9c^ - 4a^ < when a -)■ +cx). D 

Fig. Uni which is taken from Trinh's thesis shows a section of the 
surfaces Sn by the plane a = —9. Similar pictures can be seen in 

mm- 

Computational evidence suggests that each Sn has the shape of an 
infinite funnel with a sharp end stretching towards a = — oo. This 
end probably corresponds to 6 — )■ i, where b is the asymptotic value 
as in Figs. [T31IT41 A — )■ — oo as a — )■ — oo on Sn as the picture in [T3] 
suggests. For every real Oq the section of Sn by the plane a = ao is 
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Figure 16. Section of the surfaces Sn-, n = 0, 1, 2, 3 by 
the plane a = — 9. 



an oval that projects on the c-axis 2-to-l. We only proved that this 
section is compact for a large enough. For n = 0,l,2,..., the funnels 
Sn are symmetric with respect to c H- — c, S'„+i lies above 5"^ and S'^+i 
is wider than S'„,. 

Remark 4.6. In general, it is hard to say anything explicit on the cor- 
respondence between the parameters a, c in the potential and Nevan- 
linna parameter b. Some information on this correspondence can be 
extracted from symmetry and degeneration considerations. In the be- 
ginning of the proof of Theorem we noticed that the line c = corre- 
sponds to the circle \b — i/2\ = 1/2. We can determine now the sign of c 
for b inside and outside this circle. Degeneration used in the proof cor- 
responds to convergence of fe to a real non-zero point. Formula (14.1 3p 
shows that c < when e — )■ 0. So negative c correspond to the exterior 
of the circle and positive c to its interior. 

Now we state the result about the second PT-symmetric family of 
quartics. 

Theorem 4.7. The real QES part of the spectral locus of Ili4.9\ ) consists 
of \m/2\ simple disjoint analytic curves T*^^, n = 0, . . . , [771/2] — 1, 



SINGULAR PERTURBATION, PT-SYMMETRIC FAMILIES 35 

properly embedded curves which for a > project onto the ray a > 
2-to-l. When (a, A) G F*, the eigenfunction has 2n non-real zeros. 

The proof is completely similar to the proof of Theorem 14.11 

The problem of study of the whole real part of spectral locus of 
( 14. 9p . as a two-parametric family with real m seems quite interesting 
and challenging. A picture of the spectral locus for ?7i = 3 can be seen 
in g]. 



Theorem 4.8. Every equation of the form ^4-l\ ) with polynomial P 
of degree 4 which has a solution with 2n non-real zeros and infinitely 
many real zeros is equivalent by a real affine change of the independent 
variable to equation ( [^.g[ j with (a, c, A) G Sn- 

Every equation of the form Ili4.1\ ) with polynomial P of degree 4 
which has a solution with finitely many zeros is equivalent to ( [^.g| ) 
with (a, A) G T'^ „ by an affine change of the independent variable. 

The proof is similar to that of Theorem 14.41 
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